Using and Constraining Nonforward Parton Distributions 

Deeply Virtual Neutrino Scattering in Cosmic Rays 
and Light Dark Matter Searches ^ 
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Abstract 

We overview the construction of Nonforward Parton Distributions (NPD) in Deeply Virtual 
Compton Scattering (DVCS). Then we turn to the analysis of similar constructs in the weak sector 
(electroweak NPD's). We argue in favour of a possible use of electroweak DVCS (EWDVCS) as 
a rare process for the study of neutrinos in cosmic rays and for light dark matter detection in 
underground experiments. 



^ Based on a talk presented by Claudio Coriano at QCD@work 2003, Conversano, Italy, June 2003 



Nonforward Parton Distributions are an important construct in the parton model [T], [2] and appear 
in Compton Scattering (CS) (Fig. [H) in the generalized Bjorken region [H], where the process is hand- 
bag dominated. They generalize ordinary parton distributions (Fig. [2l). Proofs of factorization to all 
orders of processes of this type have also been presented |1| and there is a large amount of work (see 
for instance 0) devoted to the subject, accompanied by ongoing experimental efforts to measure the 
corresponding cross section In the electromagnetic case the measurements are difficult, given the 
presence of a dominant Bethe-Heitler background. One of the possibilities to perform the measurement 
is through the interference between the Bethe-Heitler and the DVCS hand-bag diagram by electron spin 
asymmetries. 

The analysis of these distributions, including their modeling |lZ|[IHl and the theoretical and phe- 
nomenological study of the leading-twist and higher- twist contributions in the context of ordinary DVCS 
(with a virtual or a real final state photon) 110] has also progressed steadily. In this talk we will briefly 
outline some of the features of these distributions and elaborate on the possibility to use some variants of 
them in the description of electroweak processes induced by neutral and charged currents on nucleons. 
Possible applications of these generalizations which we are suggesting are in deeply virtual neutrino 
scattering for diffuse neutrinos and in the detection of light dark matter. Results of this analysis will 
be presented elsewhere. 

2 The DVCS Domain 

A pictorial description of the process we are going to illustrate is given in Fig. El where a lepton of 
momentum / scatters off a nucleon of momentum Pi = p via a gauge boson exchange; from the final 
state a photon and a nucleon emerge, of momenta q2 and P2 == P' = p + r respectively. In the 
deeply virtual limit, analogous to the usual Bjorken limit, the final state photon is on shell, and a large 
longitudinal (light cone) momentum exchanged is needed in order to guarantee factorization. 

The regime for the study of NFPD's is characterized by a deep virtuality of the exchanged photon 
in the initial interaction {e + p ^ e + p + •y) { > 2 GeV^), with the final state photon kept on-shell, 
large energy of the hadronic system (W"^ > 6 GeV^) above the resonance region and small momentum 
transfers |t| < 1 GeV^. In the region of interest (large and small t) the Bethe-Heitler background 
(~ 1/t) is dominant and the 1/Q behaviour of the virtual Compton scattering amplitude (VCS) render 
the analysis quite complex. A dedicated study of the interference BH-VCS in order to explore the 
generalized Bjorken region is therefore required. 

2.1 Nonforward Parton Distributions 

In the case of nonforward distributions a second scaling parameter ( (^) controls the asymmetry between 
the initial and the final nucleon momentum in the deeply virtual limit of nucleon Compton scattering. 
Both the inclusive DIS region and the exclusive ERBL region can be analized with the same correlator. 
We recall that in the light-cone gauge the (off forward) distributions H{x,^) [T] is defined as 

H,{x,^,A')) = i/^e-^"^-(P'|^,(0, ^,0^)^7+^,(0, ^,Ox)|P) (1) 

with A = P' - P,A2 =t P+ = 1/2(P + P) [IJ (symmetric choice) and ^P = 1/2 A+. 

This distribution describes for x > ^ and x < the DGLAP-type region for the quark and the 
antiquark distribution respectively, and the ERBL (see (H], [H]) distribution amplitude for — ,^ < x < 
In the following we will omit the A dependence from Hq. 

The most common procedure is to use double distributions, P(x, y) defined with a symmetric choice 
of the external momenta and relate them to off-forward distributions H{v,^,t) p], incorporating in a 
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Figure 1: Compton Scattering 




Figure 2: Forward Parton Distributions 



single interval with — 1 < f < 1 both the quark and antiquark parton distributions with an asymmetry 
parameter ^, 



/-I dx' /l^l^^,, dy'6ix' + iy' - v)F{x\ y\ t) , (2) 

Below we will neglect the t dependence. The singlet/ nonsinglet decomposition of the evolution can be 
carried out as usual taking linear combination of flavours 

H^{v,0 = Ha{v,0±H3{-v,0- (3) 

The evolution can be analized in various ways. One possibility is to use nondiagonal distributions 
[H], related to the asymmetric distributions introduced by Radyushkin (see fSJ) 

\-^2 - -[^.'■^j - {i-c/2) ' y^> 

where vi G [—^,1], V2 G [—1,^]. The parameter C = f+lp+ characterizes the so-called skewedness of 
the process and is defined in the interval ( G [0, 1]. Variables ^ and ( are related in the DVCS limit by 
^ = C/(2 — C) and ( = Xbj, conciding with Bjorken variable x (modulo terms of 0{xbjt/Q^ j). 

The inverse transformations between the vs and Xs are easily found in the various allowed regions 

A'l - C/2 C/2 - A'; 

^W) = (7) 
In the gluon case one can use any either Vi or V2 equivalently. 



Figure 3: Leading hand-bag diagrams for the process 

One finds 

^-(X, C) = ^q^g^ = I\ dx' dy'6 {X' + iy' - v,) , 

^^--(X, C) = = I\ dx' dy'6 {X' + iy' - v,) . (8) 

with 1 > f 1 > —1 < f 2 < ^. The generation of initial conditions for nonforward distributions is 
usually based on factorization, as suggested originally by Radyushkin, starting from the double distribu- 
tions {x' — y' factorization). These are obtained using some profile functions 7r() which characterizes the 
spreading in y' of the momentum transfer r, combined with an ordinary forward parton distributions. 
The latter is appropriately extended to the x' < to describe antiquark components, as pointed out 
long ago ^2]- For instance, in the quark/gluon case one obtains [IB] [2] 



F^''^{x',y') = 7r''{x',y')r'^{x') 
3{l-\x'\f-y" 



4 (l-|x'|)3 
F'^ix',y') = n^{x',y')f^{x') 

15 ((1 - \x'\)^-y'y 



r^\x') 



fix') , (9) 



16 {l-\x'\y 

for quark of flavour a, where the quarks and gluon distributions are extended to x < 0, as 

f^{x) = xg{x,Qo)e{x) + \x\g{\x\,Qo)Qi-x) , 
P'\x) = g'^(x,Qo)e(x)-(g")(|a;|,Qo)e(-a;). (10) 



J^i^^{X, C) = f dx'n'^ (x', ^) g'^(x') 



(11) 



For the anti-quark, since V2 = — fi one may use eqs. (l5l8l) with t>2 — > — fi, and, exploiting the fact that 
f'^{x) = —q{\x\) for X < 0, one arrives at 



J^^XX) = dx'rc^ U',^^)r{\x'\)- 



(12) 



The expression for the gluons is similar. Performing the y' integration, one obtains in the two 
regions, DGLAP and ERBL respectively 



J'^iX, C) = llLc dx'n^ (x', yix')) x' g{x') 

^ l-C 



(13) 



and 



T^{XX) 



/(f rfxV»(x',y(x')) x'^(x') 



+ j^'^ dx'TT^ (x', y{-x')) x'g{x') 



(14) 
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Figure 4: Unitarity diagram in DIS (forward Compton scattering) 



2.2 A Comment 

Parton distributions are quasi-probabilities, very similar to Wigner functions (as first observed in (H), 
though the momentum dependence of the latter is not found easily found in the former. One possibility 
in this direction has been suggested recently [H]. The formal definition of a parton distributions is in a 
non-local light cone correlator, built through unitarity (at least in the forward case). In the nonforward 
case there is no optical theorem (Fig. Hj) that holds, and the construct is a genuine correlation function 
in the nucleon state. We have pointed out in [Zj the existence of a formal relation between evolution 
equations and kinetic equations, via a Kramers-Moyal expansion of the non-forward evolution equation. 
This relation carries a strong similarity to the well known relation between Wigner functions and their 
associated evolution equations, given by differential operators of arbitrarily high orders (Moyal products) 
as known from the phase-space approach to quantum mechanics Indeed the definition of a parton 
distribution takes place through a Wigner- Weyl transform, limited to the light cone domain. 

The evolution equations describing NFPD's are known in operatorial form [in]. Single and double 
parton distributions are obtained sandwiching the operatorial solution with 4 possible types of ini- 
tial/final states < >, < p|...|0 >, < >, corresponding, respectively, to the case of diagonal 
parton distributions, distribution amplitudes and, in the latter case, skewed and double parton distri- 
butions [U [Tl|. The ansatz for the general solution of the evolution equations up to next-to-leading 
order for Hq{x,^) is given by [Z| 

where we have introduced arbitrary scaling coefficients {Ak, Bk), now functions of two parameters (x, ^). 

The evolution can be written down in various forms. The one we have exploited comes from kinetic 
arguments on which we now briefiy elaborate. 

In the non-forward case (DGLAP-type evolution) the identification of a transition probability for the 
random walk f7] associated with the evolution of the parton distribution is obtained via the non-forward 
transition probability 



wcixly) = -Cf- 



TT y — X 



y y-C 



e{y > x) 
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Figure 5: Overlap of wave functions in proton Compton scattering 




Figure 6: ERBL factorization of proton Compton scattering 



w'Ay\x) = -Cf 



and the corresponding master equation is given by 



e{y < x) 



dr 



dywc{x\y)J^q{y,C,r) - / dyw'^{y\x)J^g{x,C,r) 

Jo 



(16) 



(17) 



that can be re-expressed in a form which is a simple generalization of the formula for the forward 
evolution 

A Kramers-Moyal expansion of the equation allows to generate a differential equation of infinite 
order with a parametric dependence on ( 



dTa 



dlogQ'' 



dyw(^{x + y\x)J^q{xXi'T)+ / dyw'Ax + y\x)J^q{xXi'T) 

a{x) JO 

- E / dy^^d/' {w^x + y\x)J^qix, C, r)) . 



(18) 



If we arrest the expansion at the first two terms {n = 1, 2) we are able to derive an approximate 
equation describing the dynamics of partons for non-diagonal transitions which is similar to constrained 
random walks [Z|. 



2.3 Possible Constraints 

Kinetic analogies are an interesting way to test positivity of the evolution for suitable (positive) bound- 
ary conditions, though their practical use appears to be modest, since it is limited to leading order 
in the strong coupling Other interesting constraints emerge also from other studies, based on 

a unitarity analysis which also may be of help. Their practical impact in the modeling of these 
distributions has not been studied at all. 
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involving iorwaru ana noniorwaru aistriuutions at tiie sanie tinie, siiouia ue anaiizeu careiuiiy in per- 
turbation theory. All the anomalous dimensions are known up to next-to-leading order and the study 
is possible. Another input comes from Ji's sum rule on angular momentum [ij, which needs to be satis- 
fied. One can view Ji's sum rule also as a constraint between the evolution of angular momentum and 
the nonforward parton distributions for a vanishing t. The evolution of all the distributions, including 
angular momentum is also known, and there is wide room to play with the boundary conditions. Other 
constraints may come from a combined sum rule analysis for proton Compton scattering, along the 
lines suggested in |Tn|. The soft mechanism based on the overlap of wave functions (see Fig. [HI), pre- 
asymptotic (Fig. Ej) may be studied using dispersion theory of 4-current correlators (Fig. EI), especially 
for scattering at wide angle. 



3 EWDVCS 

We now come to discuss an interesting generalization of the standard DVCS process. 

The electroweak version of the process (EWDVCS), in the regime that we are interested in, requires 
the notion of electroweak NPD's to be conveniently described. A charged current version of this inter- 
action, with an electron or a neutrino in the final state, the latter coming from an electron, can also 
be contemplated. We are going to argue that there are motivations for studying these extensions for 
astroparticle applications. 

If we entertain the possibility of extending DVCS to the weak sector we will be facing the issue of 
large suppression in the cross section. For charged currents there is a similar BH background, but not 
for neutral currents, for instance in neutrino-nucleon scattering, this aspect is obviously not present. 
We shoukd keep in mind that it is very likely that most of our hopes to detect dark matter using ground 
or underground based facilities, whatever they may be, have to rest considerably on our capability to 
detect weakly interacting particles of various masses and in various kinematical domains. 

Differently from standard accelerator searches, in the astroparticle case the range of these interactions 
vary from low to high to extremely high energy. 

The example of Ultra High Energy Cosmic Rays (UHECR) [2l| and the undeterminations associated 
to the estimates of their energy both in the strong sector (for instance at the Greisen-Ztsepin-Kuzmin 
(GZK) cutoff and around it for p-p first impacts) and in the weak sector - such as in the detection 
of neutrino-induced horizontal air showers - are clear examples of a theory that is unable to be fully 
predictive. We remark that arguments linked to parton saturations and to the small-x impact of the 
distributions for Ultra High Energy (UHE) neutrino scattering have been formulated only recently and 
show a considerable readjustment [21] of the UHE cross section. In particular, our overall understanding 
of this domain is still under debate, especially from the S-matrix perspective Neutrino studies have 
only been confined to the atmospheric and solar case, for obvious reasons. Though the fiuxes may 
change considerably with a varying energy, it is conceiveable that other cosmic ray experiments may be 
looking for bounds on the parameter space of the various models which have been recently proposed. 

In the search for dark matter it is essential to have an open view on the possible mass range of the 
possible candidates, including lower mass candidates j22j. While dark matter searches can be rather 
indirect and can rely on new astronomical searches with excellent perspectives for the future - such as 
weak gravitational lensing (gravitational shear measurements) or other observational methods having a 
statistical imprint - [201 weak or superweak interactions are still likely to be studied extensively using 
nucleon targets. We believe that the few GeV region of QCD has very good chances to be relevant in 
the astroparticle territory, and for a just cause. One additional photon in the final state plus a recoiling 
nucleon is a good signal for detecting light dark matter and both the fiuxes and the cross sections should 
be accurately quantified [26] . 
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Figure 7: Dispersive Approach to proton Compton scattering using current correlators (above). The 
spectral density of pion Compton scattering (below). 

4 Conclusions 

We have illustrated some of the basic features of NFPD's and we have emphasized that there are other 
possible applications of the theory, especially in the weak sector, which may be relevant in the astropar- 
ticle domain. Any attempt to obtain a complete coverage of neutrinos and dark matter interaction with 
nucleons in underground experiments should also include searches in the intermediate energy region of 
QCD with far reaching consequences. 
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